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Abstract. - Recently we have proposed a s + g-wave model for the superconductivity in 
borocarbides, YNi2B2C and LuNi2B2C 0IHini- In the present paper we first summarize ther- 
modynamic properties of s + g-wave model. Then we shall analyse the recent Raman spectra 
data of RNi2B2C (with R=Lu and Y) by Yang et al. [J. The present model appears to describe 
salient features of the Raman spectra. 



Introduction. - The Superconductivity in rare earth borocarbides is of great interest 
[Slini . In particular its interplay with magnetism and superconductivity is very fascinating . 
However in the following we limit ourselves to superconducting borocarbides LuNi2B2C and 
YNi2B2C. They have relatively high superconducting transition temperature Tc = 16.5 K 
and 15.5 K, respectively. Although the dominance of s-wave component in A(k) has been 
established by substituting Ni by Pt and subsequent opening of the energy gap the 
superconductivity exhibits a number of peculiarities common to nodal superconductivities 
|1(JI11| . For example both the ^/H dependence of the specific heat and the H linear dependence 
of the thermal conductivity indicate that the superconductivity has the nodal excitations 
|12lll3ll4lll5iri 6 17, 18 similar to d-wave superconductors in high Tc cuprate superconductors. 
Further the presence of de Haas van Alphen oscillation in the vortex state of LuNi2B2C 
down to H — 0.2i7c2 indicates again the nodal superconductors ^Hl- In conventional s-wave 
superconductor de Haas van Alphen oscillation would disappear for H < O.SHc ,20,. Further 
the upper critical field determined for LuNi2B2C and YNi2B2C [2] in a magnetic field within 
the a-b plane exhibits clear fourfold symmetry reminiscent to d-wave superconductors |22| . 

These experiments indicate clearly A(k) in borocarbides has to have an anisotropic s-wave 
order parameter. Further a) A(k) has to have the nodal structure or the quasiparticle density 
of states, N{E) ^ \E\ for \E\ <C A where A is the superconducting order parameter (i.e. 
the maximum of A(k)), which gives both the v^dependence of the specific heat and the H 
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Fig. 1 - The specific heats normahzed by the normal state one are shown for s + p-wave, s-wave and 
d-wave superconductors. 7jv is the Sommerfeld Constant. 



Fig. 2 - The spin susceptibility normalized by the normal state are shown for s -I- (j-wave, s-wave and 
d-wave superconductors. 



linear thermal conductivity in the vortex state, b) the nodal structure has to have the fourfold 
symmetry within the a-b plane and to be consistent with the tetragonal symmetry. These two 
constraints appear to suggest almost uniquely 

A(k) = ^A(l -sin4 6'cos(4<^)) (1) 

or s+g-waye superconductor. Here 9 and </> are the polar and azimuthal angles describing k. 
Contrary to Ref. we have minus(-) sign in front of the g-wave term. This corresponds to 
point nodes at [100], [010], etc. Those positions of the nodal points are consistent with the 
magnetothermal conductivity data More recently the point node at the same positions 
have been seen by the magnetospecific heat measurement in the vortex state of YNi2B2C by 
Park et al. We shall see later the Raman spectra data from both YNi2B2C and LuNi2B2C 
are consistent with Eq.(l). Recently the above Raman data have been analysed in terms of 2D 
s + g-wa.ve model by Lee and Choi [21] • Unlike the present model their model has line nodes. 
Therefore their model cannot describe the magnetothermal conductivity data 2] . Further the 
description of the Raman spectra within this model is unsatisfactory. 
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Thermodynamics. - First we determine the temperature dependent energy gap A(r) 
within the weak coupUng theory. Here A(r) is the maximum value of the energy gap A(k). 

f / dE (Re , ^ = > tanh(^) (2) 

where / = ^(l — sin^ 6* cos(4(/))). A and Ec are the dimensionless coupHng constant and the 
cut-off energy(:» A), respectively. Further < . . . >= / dQ/in. 
In the vicinity of T = Tc, Eq.(2) gives 



2 2(27rT)2</2> T 

7C(3) <r> ^-'''Yj 

with Tc = ^EcC^^/^ and 7 = 1.78107 ... the Euler constant. On the other hand for T < Aq 



where 



= 1.9635(|)3 + 14.2(J)^.. (5) 



Ao = A(0) = ^T^exph < f >-^< fln\f\ >] ~ 2.76Te (6) 

TT 



Then the specific heat Cs is given by 



.Hz!„„T{|c(3)(I) + g(I)^+...) (8) 



where E = y^S,"^ + A^/^ and A^o is the density of states in the normal state. In Fig.l we 
show Cs/{^^NqT) versus T/Tc for s + g-wave, s-wave and d-wave superconductors. As is 
readily seen the specific heat of s + g-wave superconductor is very similar to the one in d-wave 
superconductors . Also the spin susceptibility of s -f g-wave superconductor is shown in 
Fig. 2, which is very similar to the one in d-wave superconductors. 

The superfluid density, on the other hand, has the axial symmetry. The superfluid density 
in the a-b plane is given by 

"-^"^^^ - l-± rdC< sin^ sech^d;) > (9) 



^(ln2)(I)~— ... (10) 



while the superfluid density parallel to the c-axis is given by 



o poo rp 

l-^y^ d^<cos^O sech\ — )> (11) 
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Fig. 3 Fig. 4 

Fig. 3 - the superfluid density in the a-b plane for s + g- and d-wave superconductors. 

Fig. 4 - the superfluid density parallel to the c-axis for for s + g- and d-wave superconductors. 



In Fig. 3 and Fig. 4 we show the superfluid density in the a-b plane and parallel to the c- 
axis, respectively. The superfluid density in the a-b plane is very similar to the one in d-wave 
superconductor. Further the superfluid density parallel to the c-axis is somewhat similar to 
the one in the d-wave superconductor which is due the coherent Josephson tunneling j26| . We 
note also T^^ ~ behavior has been already observed in Ref. j27| . 

Further the thermal conductivity within the a-b plane exhibits the universal heat conduc- 
tion |2HlEni 

^/T (13) 

' 8 Aom ^ ^ 

in the limit K. Also the thermal conductivity for H || c gives 

^(H)/.. = (14) 

for T, VFA <C VaVeH, where F and Va are the electron scattering rate and the Fermi velocity 
within the a-b plane, respectively. k„ is the one in the normal state and Va is the Fermi 
velocity within the a-b plane. Indeed the _ff-linear thermal conductivity is observed recently 
by Boaknin et al. |15j . 

Also very recently the specific heat of the vortex state in YNi2B2C in a magnetic field 
within the a-b plane is observed by Park et al. 23 . It exhibits cusps at H || a and H || b 
typical to the point nodes in s + g-wave superconductors. 

Electronic Raman Spectra. ~ We consider the case the polarization vector of photons lie 
in the a-h plane. Then the Raman spectra are given by |3()ll31j 

5,(a;/2A) = Im[ < 7,f A > - ^^y^] (15) 

where 7Aig = A/2cos(40fc), 7^19 = \/2cos(20fc) and 7529 = V2sin(20fc) and 

A = A' + a" (16) 
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Fig. 5 - The theoretical Raman spectra in Aig, Big and B2g modes using s + g model at T = K 
(left panel) and the experimental Raman spectra for YNi2B2C taken at T = 6 K (right panel) are 
shown. 



A' = 



2 



tan"^( 



X" = 



2x - p 



(17) 
(18) 



where x = u}/2A and / = A(k)/A and < . . . > means the k(angle) average for the Fermi 
surface. 9{x) is the step function, i.e., 9{x) = 1 for x > and 9{x) = for x < 0. We note for 
the present model the second term in Eq.(lO) docs not contribute for Big and B2g modes due 
to the symmetry constraints. The electronic Raman spectra from these 3 modes are obtained 
numerically and shown in Fig. 3 (the left-side panel). In parallel to the theoretical result we 
show the experimental data for YNi2B2C taken at T = 6 K. 

We note Aig mode is very consistent with the observed spectra. Also the low frequency 
parts {to < A(T)) of both the Big mode and the B2g mode are very consistent. On the other 
hand the theoretical curve for the B2g mode exhibits a cusp at uj = 2A{T). whic;li is not 
seen experimentally. We don't know if the cusp-like feature will disappear with increasing 
temperature or not. Also the peak position of the Big mode is somewhat in the higher energy 
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than that of the B2g mode. Again we don't know if this is the effect of the temperature. In 
any event we may conclude that s + g-wave model captures the main feature of the Raman 
spectra of YNi2B2C. Also the peak position in the B2g mode give A(T). 

The weak coupling theory gives Ao for YNi2B2C(Tc ~ 15.3 K) and LuNi2B2C (T^ ~ 15.7 
K) is 42.2 K(3.64 meV) and 43.3 K(3.73 meV), respectively. On the other hand the data 
at T/Tc ~ 1/2 indicate Ao=50.4 K and 64.7 K for YNi2B2C and LuNi2B2C, respectively. 
Therefore we may conclude that the borocarbides superconductors are in the intermediate 
coupling region. However, clearly a further experiment at low temperatures is highly desirable. 

Summary. - We have analysed further the s + g-wave superconductors proposed in [2|2] . 
This model appears to describe the recent specific heat data as well. We have also studied 
the Raman spectra reported in Ref. jl] . The present model appears to capture the main feature 
of the observed spectra. This further confirms the presence of the order parameter A(k) with 
point nodes at k =(100), (010), (100), and (010). 
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